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Abstract. We review some major, open issues in the current modelling of low and intermediate 
mass, main sequence stars based on seismological studies. In the present paper, the solar case is 
discussed together with current problems that are common to the Sun and stars with a structure 
similar to that of the Sun. Several additional issues specific to main sequence stars other than 
the Sun are reviewed and illustrated with a few stars observed with CoRoT in a companion 
paper. 



1. Introduction 

After more than two decades of helioseismology, almost four years of asteroseismology with 
CoRoT [1] and almost two years of intensive asteroseismology with Kepler [2], we review some 
major, current open issues about the internal structure of the Sun and solar-like oscillating 
stars. We discuss here the solar case, this also applies to oscillating stars that have a similar 
internal structure. For sake of brevity, we decided to include only unpublished figures and to 
cite published figures in the text. Several recent reviews exist on the topic, for instance [3], [4], 
[5]. 



2. The Sun 

As it is well known, the Sun is a particular case. It is the closest star, and as a result we 
know with a high precision the luminosity, mass (through the product GMq), radius, age and 
individual surface abundances of chemical elements^. Furthermore, a wealth of very accurate 
seismic constraints are available and have been successfully used. Inversion of a large set of mode 
frequencies has provided crucial information on the structure of the Sun, see for instance [5]; 
[6]. Accordingly, the following constraints must all be satisfied by any calibrated solar model: 
radius at the base of the convective envelope rf,zc, surface helium abundance Yg, sound speed 
profile c(r), internal rotation profile and location of ionization regions. To some extent, these 
constraints are found to be independent of the reference model [7]. The current major challenges 
and open issues in the solar case then are: 

• what are the values of the surface abundances, more specifically the oxygen abundance? 



^ although some of these latter are still debated, see below 



Table 1. Value of the solar photospheric metallicity from 1993 to 2010 from Table 4 in Asplund 
et al [10] complemented by Caffau et al [11] 

GN93 GS98 AGS05 AGS09 Lod09 CafflO 
Z/X 0.0244 0.0231 0.0165 0.0181 0.0191 0.0209 



• what is the origin of the discrepancy between the seismic sound speed and that given by 
models below the convection zone? 

• what are the dominant physical mechanisms responsible for a uniform rotation in the 
radiative region and a differential rotation in the convective zone? It is worth noting that 
this is the opposite in current calibrated ID solar models: the convection zone is assumed to 
rotate uniformly and the rotation profile in the radiative zone is found to vary with radius 
if waves and/or magnetic field are not taken into account (see below) ! 

• how can we model properly near surface layers and the convection-pulsation interaction? 

• how to succeed in probing the core? 

• how to model oscillation mode line widths and amplitudes? 

Some of these uncertainties about the Sun have consequences on the modelling of stars other 
than the Sun. On the other hand, some problems that are encountered with seismological studies 
of stars other than the Sun can be studied first with our well known Sun. These two points of 
view are discussed below. 

2.1. Initial abundances: the solar mixture 

In the late nineties, with the GN93 solar abundances [8], the seismic sun and calibrated solar 
models were in agreement by 1 to 5% for the sound speed profile and location of the base of 
the convective zone [9]. However between 1993 and 2010 several revisions of the photospheric 
solar mixture were performed. Noteworthy, 3D model atmospheres including NLTE effects and 
improved atomic data have led to a substantial decrease of the C, N, O, Ne, Ar abundances and 
in turn to an important decrease of the solar metallicity [Z/X) (see Table 1). 

Today, an internal consistency of the abundance determination from different ionization 
levels of a given element seems to have been obtained and a consensus between independent 
determinations seems to be reached (i.e. this shows the unvaluable benefit of independent 
approaches). However the revised modern solar chemical composition [10] leads to a strong 
disagreement between the sound speed of a calibrated solar model and that of the seismic Sun. 
This is also true for the depth of the convection zone and the envelope helium abundance (see [5] , 
for a review). Fig.l shows the difference of the sound speed profile between Basu's seismic model 
[12] and several models computed with Cesam2k [13] using some of the abundance mixtures listed 
in Tab.l. This figure shows that the discrepancy at the base of the convection zone is strongly 
increased when the revised abundance mixture AGS05 is used. This is mainly due to the decrease 
in oxygen and neon and in {Z/X)s which consequence is a decrease of the radiative opacities. 
With the newly derived abundance mixture AGS09 which has achieved some consensus, the 
discrepancy slightly decreases but remains significant. 

2.1.1. Possible origins of the discrepancy First there could be some errors in the opacity 
derivation. As a check, the comparison between OPAL and OP opacities for a calibrated solar 
model shows that OP opacities give a (slightly) better fit than OPAL. However a change in 
opacity of about 30% at the base of the convection zone, at a temperature of 2 10^ K would be 
necessary to compensate for the effects of the change in mixture [15, ?] and according to [16] 
there appears to be no way to change the OP opacity by such an amount (see also [17] and [18]). 




Figure 1. Differences in the sound speed profiles for the Sun between Basu's seismic model and 
calibrated solar models using various abundance mixtures and different assumptions for the physical 
description. All models have been calculated with the Cesam2k code except model S [14] which is based 
on the old solar mixture. 



Other possible causes have been discussed but none of the related improvements can reconcile 
simultaneously all the seismic constraints listed above [19], [20], [21]. 

In order to reconcile the seismic and theoretical sound speed, one needs higher opacities or 
higher helium below the upper convective zone (UZC) i.e. a higher He gradient. Any mixing 
below the UZC which smooths the gradient therefore would go in the wrong direction. In that 
case, one would rather need an advection process. It is not clear whether internal waves and/or 
which hydrodynamical instabihty could act in the sense of increasing helium below the UZC. 
Other further developements are underway such as taking magnetic effects into account in the 
derivation of solar abundance corrections. An investigation based on 3D magnetoconvection 
simulations shows that corrections to the solar abundance can be significant [22]. 

Abundances of other stars are determined by reference to the Sun, hence all stars are affected 
by errors or by inaccuracies in the solar mixture. Can other stars be discriminating in that issue? 
Note that one must wonder about the impact of the inconsistency which arises when modeling 
other stars using solar mixtures inferred from 3D model atmospheres if the stellar [Fe/H] itself 
has not been determined from a dedicated 3D model atmosphere. 



2.2. Nuclear reaction rates 

In stars, (most commonly charged induced) nuclear reactions occur at low energies (10- 300 
keV) in the Gamow peak -which corresponds to the maximum probability of the reaction. The 
cross sections that are governed by Coulomb barriers and resonances show a strong and complex 
dependance on energy and globally decrease steeply towards low energy. 
Cross sections are usually written in the form 

aiE) = ^ e-'-^ (1) 

where S{E) -the astrophysical factor- contains everything that concerns the nuclei and nuclear 
physics and varies slowly with E. The exponential term is related to the Coulomb barrier and 
relative velocity of the nuclei. Experiments which provide measurements of the cross sections in 
the laboratory generally occur at higher energy than the Gamow peak. Extrapolation of S{E) 



to lower E is then necessary. This is difficult to achieve and cannot take into account possible 
unknown resonances occuring at low E. When the S-factor cannot be measured, the reaction 
rates are obtained from pure theory. Recent significant progress in laboratory and theory (hence 
better determination of S-factors down to the Gamow peak) has been achieved as discussed in 
the comprehensive reviews by [23], [24], see also [25]. 

2.2.1. Hydrogen burning reaction rates: Uncertainties still exist for the pp chain and CNO 
cycle reaction cross sections in the Sun. They are due to the difficulty to estimate the S-factor at 
low energy and to determine the exact role of the electron screening (both in the laboratory and 
in the star). For most of the concerned reactions, S"- factors are determined from extrapolation 
of experimental data to low E but now, for some key reactions, energies corresponding to the 
Gamow peak are accessible through measurements provided by the LUNA experiment at Gran 
Sasso [26]. This is the case of the ^He(^He,2p)'^He reaction (see Fig. 1 in [26]) for which 
the 5'(£^)-value has not been significantly modified but the error bars are largely decreased. 
Furthermore the cross section of the N^^{p,^)0'^^ reaction -the leading reaction in the CNO 
cycle- has now been measured down to energies relevant for IMq stars on the red giant branch 
(see Fig. 19 in [26]) and a resonance is observed. This represents a significant advance although 
an extrapolation is still needed towards solar conditions. Noteworthy, the revision is important 
with a decrease of S{E) by 50%. This has crucial consequences for solar and stellar structure. 
In the solar core, the CNO cycle efficiency is reduced from 1.6 to 0.8% of the total energy. For 
main sequence stars slightly more massive than the Sun the occurrence of a convective core 
depends on the energy production. With the new Af^^(p, 7)0^^ rate the convective core appears 
at higher stellar mass and is less massive for a given mass (see for instance the case of a 1.2 
Mq low metallicity star in Fig. 14 by [27]). Finally cross sections that are obtained from pure 
theory can be constrained by helioseismology [28], [29], [23]. 

2.2.2. Electron Screening Electron screening is based on the Salpeter's formula [30] with the 
underlying physical picture: the cloud of electrons decrease the repulsive Coulomb effect between 
interacting nuclei with the result of a decrease of the Coulomb potential and an enhancement 
of the reaction rate. This is a static description. It is currently not clear whether dynamic 
effects of the interacting ions can significantly change the impact of the screening in nuclear 
interactions for stars and wherther it must therefore be taken into account. The energy that 
initially fast moving interacting ions have when they get close enough to interact can be lower 
than the mean value of the medium; accordingly their reaction rate is reduced compared to 
Salpeter's prescription. This effect is difficult to quantify and relies on results from numerical 
simulations [31], [32], [33], [34] and the issue is not settled yet. [29] looked at the impact of 
changing the electron screening compared with the classical Salpeter's formula on the solar 
sound speed profile. They found that the solar seismic constraints do not allow variations larger 
than 1% when using GS98. [3] computed a solar model equivalent to the S model but assuming 
the extreme case of no screening at all. Their Fig. 12 shows the difference in the sound speed 
profile between the seismic sun [12] and model S but switching off e- screening. A decrease 
of the reaction rate by swithcing off the e-screening clearly increases the discrepancy between 
observation and model. This is in agreement with previous results of [23] who showed that the 
large sound speed discrepancy with the use of AGS05 mixture is reduced to the previous level 
obtained with GS98 when the pp-reaction and the screening factor are increased up to 15% 
(Fig. 3 in [23]) In that case, the surface helium and the depth of the convective envelope are 
in agreement with seismic determinations but the sound speed profile in the core significantly 
deviates from the seismic solar sound speed. 



2. 3. Transport of chemicals and angular momentum 

Rotationally induced transport The physical origin of the uniform rotation profile in the radiative 
zone of the Sun unravelled by helioseismology is still debated. Rotationally induced transport 
of angular momentum resulting from a competition between shear-induced turbulence and 
meridional circulation driven by surface angular momentum losses is not able to make the 
rotation uniform in the radiative region of the Sun. For details, the reader is refered to [35]. The 
impact of rotationally induced transport on the solar sound speed profile has been investigated 
by [36]; [37]. [38] have computed a solar model including rotationally induced transport with 
two diff'erent assumptions about the initial velocity (slow or 'fast' sun). Their conclusion is that 
for an initially slow rotation, the microscopic difi'usion dominates whereas for an initially rapid 
enough rotation, meridional circulation dominates over turbulent shear. However in both cases, 
the discrepancy below the UCZ increases (see Fig. 9 of [38]). This is also illustrated in Fig.l 
which shows the differences in sound speed profile for calibrated solar models computed with 
AGS05 mixture computed assuming either no rotationally induced transport, or rotationally 
induced transport with no surface angular momentum loss or assuming rotationally induced 
transport with surface angular momentum loss. Adding rotationally induced transport makes 
the discrepancy with the observed sound speed worse. Indeed this process smoothes the helium 
gradient below the UZC whereas the observations seem to require an increase of helium in the 
radiative zone below the UCZ. Several prescriptions have been derived for the turbulent transport 
coefficients involved in rotationally induced transport but validitation of these prescriptions 
remain to be done [39]. 

Internal wave induced transport in radiative zones must exist in stars as these waves are generated 
at the interface between convective and radiative regions. They have been shown to transport 
angular momentum efficiently enough to make the rotation of the Sun in its radiative part rigid 
[40] . One the main uncertainties is related to the absence of a viable quantitative description of 
the generation of waves. 

Magnetic induced transport Instabilities driven by the interaction between rotation and a 
magnetic field could be responsible for the transport of angular momentum and the rigid solar 
rotation [41]. This was shown by [37] in the solar case. 

Whether all 3 processes work together to shape the solar rotation profile in the radiative zone 
or only one or two are dominant is not settled yet. 

Investigations of the impact of rotationally and wave induced transport on the structure of 
stars other than the Sun have been performed for instance to explain the Li dip [42]; [43]. [44] 
studied the impact of including both rotationally induced transport and magnetic field on the 
structure of solar like stars while [45] estimate the seismic consequences and find that with 
the type of dynamo they assume in the radiative zone, the efficiency of rotational mixing in a 
radiative zone is significantly decreased and seismic parameters are then similar to those of a 
non-rotating, non-magnetic star. 

2.4- Near surface layers 

A direct comparison of the observed and numerical frequencies of the Sun shows systematic 
differences that remain small at low frequency but increase with increasing frequency [46]. 
Several causes contribute to this offset with more or less importance [47]. They are collectively 
referred to as near-surface effects (for reviews, see for instance [48], [4]). 

2.4-.1. Surface turbulent convection One important contribution to the differences between 
the observed and calculated frequencies comes from current modelling of the outer turbulent 
convective layers of the Sun. The description of the convective outer layers of the Sun in ID 
stellar models remains quite approximate due to our inability to represent and implement in a 
ID code a 3D multiscale nonlocal process such as turbulent dynamics [49]. At the solar surface, 
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Figure 2. Echelle diagrams computed with GOLF solar data (blue dots) and using a calibrated solar 
model computed with either an Eddington (left panel) or a Kurucz (middle panel) T(t) law (red dots) 
. Filled symbols correspond to radial modes. Right panel: An additional echelle diagram obtained 
from frequencies of the solar model (red dots) but corrected for near surface effects (green dots) agrees 
with the observations (blue dots) over the fitted interval. 

turbulent convection is inefficient and therefore strongly dependent on the free parameters 
entering the local, ID formulation for convection as well as many other assumptions in the 
formulation. A comparison between frequencies computed with two of the available formulations 
(MLT and CM [50]) for instance shows that the frequency differences increase with frequency 
and reach up to 0.15% at a frequency z/ = 4. mHz i.e. a frequency offset of about 6 — 10/iHz 
([9])- 

Patched models, that is ID stellar models where the outer layers have been replaced in a 
proper way by those obtained with 3D numerical simulations, lead to frequencies in much better 
agreement with the observations [51], [52], [53], [54]. For reviews, see [55], [56]. This approach is 
valuable for studies of individual stars provided great care is taken in the patching procedure, but 
it cannot be used in a systematic investigation of a large number of stars. Indeed 3D numerical 
simulations with the required quality are quite numerical time consuming and therefore not 
available in the whole range of effective temperature, gravity and chemical composition. 

Atmosphere as boundary condition: Another longstanding problem is the description of the 
atmosphere as boundary condition when computing a stellar model [57]. In the solar case, 
one can use an empirical atmosphere derived from observations (HSRA) [58] which represents 
accurately enough the Sun atmospheric properties. But in other stars, one must rely on 
atmospheres built using a temperature-optical depth, T(r), law. Commonly used T{t) laws are 
the Eddington law or more realistic Kurucz model atmospheres [59]. In order to illustrate the 
impact of the atmosphere boundary condition, Fig. 2 (left panel) shows an echelle diagram built 
from observed solar oscillation frequencies based on SOHO data. An echelle diagram computed 
for a calibrated solar model where the temperature stratification of the atmosphere is assumed 
to follow a Eddington T{t) law is also shown (middle panel). Largest differences are seen at 
high frequency. They are significantly larger than the observational errors. The discrepancy at 
high frequency decreases (roughly by half) when a more realistic Kurucz model is used (middle 
panel) . However one needs extensive grids of such models along an evolutionary track and for 
different masses. Furthermore, model atmospheres suffer from physical imperfections and are not 
applicable over the entire range of needed masses and ages (i.e. gravity, effective temperature) 
[60]. 



Nonadiabatic effects: Another source of uncertainty comes from what is usuahy referred to as 
nonadiabatic effects. They include the effects of interaction of the wave with the radiation and 
with the convection (see for instance [48]). 

All these imperfections concur to generate significant errors on the computed oscillation 
frequencies which in the solar case amount up to 10 /xHz at high frequency. As near surface 
effects cannot yet be reliably included in theoretical frequency computations, an alternative 
has been proposed which consists in removing these effects from observed frequencies. In the 
solar case, this has been quantitatively assessed with the comparison of theoretical and observed 
frequencies of modes with small to large degrees. 

2.4-2. Correcting frequencies for near surface effects In the solar case, due to the large number 
of different modes, near surface effects can be removed from the frequencies. However, other 
procedures must be found for stars other than the Sun. Several approaches have been proposed. 
One consists in comparing the ratio of frequency combinations rather than frequencies themselves 
for low degree modes. In that case indeed near surface effects are at least partially cancelled 
[61]. On the other hand, in order to be able to use absolute frequencies, [62] proposed a means 
to correct observed frequencies for near surface effects. They first showed that the systematic 
offset between the observed and theoretically computed frequencies of the Sun is well fitted with 
a power law 

Vm.od{n) = Vobs{n) - a{Uohs{n)/h'maxf (2) 

with a, h fit to the data and Vmax = 3100 Hz is a reference frequency that corresponds to 
the frequency of maximum power in a power spectrum. The frequencies fof,s(ra) and I'rnodi'^) 
respectively represent the observed and model frequencies of radial modes with radial order n. 
For the Sun, such a correction leads for instance to small separations (io2 which yield a solar age 
consistent with the meteoritic age [63]. 

Fig. 2 (right panel) illustrates the effect of correcting / = frequencies from near surface effects 
according to Eq.2 with r = 1, a = —3.0, b = 4.9 for the same solar model (with Kurucz T(t) law) 
as in the middle panel. The corrected frequency echelle diagram coincides with the observations 
over the frequency interval that was used to derive the parameter values a, b. The next question 
is: how much the parameters a,b, vq do depend on the adopted (Kurucz or another) model 
atmosphere? Of course one must also keep in mind that the values of the fitted parameters are 
valid only over the fitted observed frequency domain (see Fig. 2 left). 

Whether the above procedure can be applied to stars that are different from the Sun is an 
open question. Elements of the answer can be obtained with the use of available 3D numerical 
simulations and ID patched models. Frequency differences for radial modes between a solar 
patched model and a non patched model are displayed in Fig. 3 (left) as a function of the radial 
mode frequency of the patched model scaled to the frequency fmax,© given according to [64] 
scaling law. The curve a{vpatch / T^max)^ is also plotted and coincides with the observations over 
the frequency interval that was used to derive the a, b parameter values over an interval up to 
^patch/^max ~ 1-2. Similar plots are displayed in Fig.3 (right) for 3 other models with different 
effective temperature, gravity or chemical abundance. For these models, we take again Umax 
according to [64] and T^ff^Q and Vmax,G from our solar patched model. The parameters a, b are 
adapted to fit the frequency differences over the largest Vpatch / i^max interval. Again a power law 
can fit the frequency differences over an interval up to I'patch/^max 

f» 1.2 - 1.4. We note that 

the frequency differences for the 3 models behave differently at high frequency than the solar 
ones. This is likely due to their lower gravity, higher effective temperature. Indeed, for these 
models, the ratio of turbulent to total pressure is higher than for the solar model (the higher 
the effective temperature, the lower the gravity, the higher the ratio Pturb/Ptot) and this results 
in larger differences between patched and non patched models. Note that for the three models 




Figure 3. left: Frequency differences for radial modes Vpatch ~ i^nopatch between a patched and a non 
patched solar models (blue dots) as a function of Vpatch/ ^max- The black dashed line represents the 
power law according to Eq. 2 with Vmax = 3.207mHz,r = l,a = — 15,& = 4.9 (For the solar models: 
Teff = 5780, g = 4.46). right: Same as left for stellar models with different effective temperature, 
gravity as labelled in the left corner. Magenta and red models have a solar metallicity whereas the 
cyan models have a lower metallicity [Fe/H] = —1. Coloured dashed lines represent the power law 
according to Eq. 2 for each corresponding frequency difference. Parameter values are Vmax = 1.462mHz, 
r = 1, a = -17, h = 3.9 (red); v,na.x = 1.098mHz, r = 1, a = -10, 6 = 4.9 (violet); and v,nax = 1.071mHz, 
r = 1, a = —15, b = 4.5 (cyan) 



(all three hotter than the solar model), the frequency differences show an oscillating behavior 
in function of the scaled frequency. Whether this behavior is real or artificially introduced by 
the patching process is not known yet. In any case, at lower frequency, the mean variation with 
frequency is well reproduced by a power law up to Vpatch / I'max ~ 1-2 — 1.4 

3. Conclusion: 

Despite the huge amount of information provided by helioseismology about the internal structure 
of the Sun, several important open issues remain that we have reviewed. They can also impact 
our understanding and modelling of stars that have a similar structure to that of the Sun. 
An increasing set of such stars are being observed by CoRoT and Kepler, and seismic tools 
developped for studying the Sun are now being adapted to study other stars. 
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Abstract. We review some major, open issues in the current modelling of low and intermediate 
mass, main sequence stars based on seismological studies. In the present paper, the solar case is 
discussed together with current problems that are common to the Sun and stars with a structure 
similar to that of the Sun. Several additional issues specific to main sequence stars other than 
the Sun are reviewed and illustrated with a few stars observed with CoRoT in a companion 
paper. 



1. Introduction 

After more than two decades of helioseismology, almost four years of asteroseismology with 
CoRoT (1) and almost two years of intensive asteroseismology with Kepler (2), we review some 
major, current open issues about the internal structure of the Sun and solar-like oscillating 
stars. We discuss here the solar case, this also applies to oscillating stars that have a similar 
internal structure. For sake of brevity, we decided to include only unpublished figures and to 
cite published figures in the text. Several recent reviews exist on the topic, for instance (3), (4), 
(5). 

2. The Sun 

As it is well known, the Sun is a particular case. It is the closest star, and as a result we 
know with a high precision the luminosity, mass (through the product GMq), radius, age and 
individual surface abundances of chemical elements^. Furthermore, a wealth of very accurate 
seismic constraints are available and have been successfully used. Inversion of a large set of 
mode frequencies has provided crucial information on the structure of the Sun, see for instance 
(5); (6). Accordingly, the following constraints must all be satisfied by any calibrated solar 
model: radius at the base of the convective envelope r^^c; surface helium abundance Yg, sound 
speed profile c(r), internal rotation profile and location of ionization regions. To some extent, 
these constraints are found to be independent of the reference model (7). The current major 
challenges and open issues in the solar case then are: 

• what are the values of the surface abundances, more specifically the oxygen abundance? 
^ although some of these latter are still debated, see below 



Table 1. Value of the solar photospheric metallicity from 1993 to 2010 from Table 4 in Asplund 
et al (10) complemented by Caffau et al (11) 

GN93 GS98 AGS05 AGS09 Lod09 CafflO 
Z/X 0.0244 0.0231 0.0165 0.0181 0.0191 0.0209 



• what is the origin of the discrepancy between the seismic sound speed and that given by 
models below the convection zone? 

• what are the dominant physical mechanisms responsible for a uniform rotation in the 
radiative region and a differential rotation in the convective zone? It is worth noting that 
this is the opposite in current calibrated ID solar models: the convection zone is assumed to 
rotate uniformly and the rotation profile in the radiative zone is found to vary with radius 
if waves and/or magnetic field are not taken into account (see below) ! 

• how can we model properly near surface layers and the convection-pulsation interaction? 

• how to succeed in probing the core? 

• how to model oscillation mode line widths and amplitudes? 

Some of these uncertainties about the Sun have consequences on the modelling of stars other 
than the Sun. On the other hand, some problems that are encountered with seismological studies 
of stars other than the Sun can be studied first with our well known Sun. These two points of 
view are discussed below. 

2.1. Initial abundances: the solar mixture 

In the late nineties, with the GN93 solar abundances (8), the seismic sun and calibrated solar 
models were in agreement by 1 to 5% for the sound speed profile and location of the base of 
the convective zone (9). However between 1993 and 2010 several revisions of the photospheric 
solar mixture were performed. Noteworthy, 3D model atmospheres including NLTE effects and 
improved atomic data have led to a substantial decrease of the C, N, O, Ne, Ar abundances and 
in turn to an important decrease of the solar metallicity {Z/X) (see Table 1). 

Today, an internal consistency of the abundance determination from different ionization 
levels of a given element seems to have been obtained and a consensus between independent 
determinations seems to be reached (i.e. this shows the unvaluable benefit of independent 
approaches). However the revised modern solar chemical composition (10) leads to a strong 
disagreement between the sound speed of a calibrated solar model and that of the seismic Sun. 
This is also true for the depth of the convection zone and the envelope helium abundance (see 
(5), for a review). Fig.l shows the difference of the sound speed profile between Basu's seismic 
model (12) and several models computed with Cesam2k (13) using some of the abundance 
mixtures listed in Tab.l. This figure shows that the discrepancy at the base of the convection 
zone is strongly increased when the revised abundance mixture AGS05 is used. This is mainly 
due to the decrease in oxygen and neon and in {Z/X)s which consequence is a decrease of the 
radiative opacities. With the newly derived abundance mixture AGS09 which has achieved some 
consensus, the discrepancy slightly decreases but remains significant. 

2.1.1. Possible origins of the discrepancy First there could be some errors in the opacity 
derivation. As a check, the comparison between OPAL and OP opacities for a calibrated solar 
model shows that OP opacities give a (slightly) better fit than OPAL. However a change in 
opacity of about 30% at the base of the convection zone, at a temperature of 2 10^ K would 
be necessary to compensate for the effects of the change in mixture (15; 3) and according to 
(16) there appears to be no way to change the OP opacity by such an amount (see also (17) 
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Figure 1. Differences in the sound speed profiles for the Sun between Basu's seismic model and 
calibrated solar models using various abundance mixtures and different assumptions for the physical 
description. All models have been calculated with the Cesam2k code except model S (14) which is based 
on the old solar mixture. 



and (18)). Other possible causes have been discussed but none of the related improvements can 
reconcile simultaneously all the seismic constraints listed above (19), (20), (21). 

In order to reconcile the seismic and theoretical sound speed, one needs higher opacities or 
higher helium below the upper convective zone (UZC) i.e. a higher He gradient. Any mixing 
below the UZC which smooths the gradient therefore would go in the wrong direction. In that 
case, one would rather need an advection process. It is not clear whether internal waves and/or 
which hydro dynamical instability could act in the sense of increasing helium below the UZC. 
Other further developements are underway such as taking magnetic effects into account in the 
derivation of solar abundance corrections. An investigation based on 3D magnetoconvection 
simulations shows that corrections to the solar abundance can be significant (22). 

Abundances of other stars are determined by reference to the Sun, hence all stars are affected 
by errors or by inaccuracies in the solar mixture. Can other stars be discriminating in that issue? 
Note that one must wonder about the impact of the inconsistency which arises when modeling 
other stars using solar mixtures inferred from 3D model atmospheres if the stellar [Fe/H] itself 
has not been determined from a dedicated 3D model atmosphere. 



2.2. Nuclear reaction rates 

In stars, (most commonly charged induced) nuclear reactions occur at low energies (10- 300 
keV) in the Gamow peak -which corresponds to the maximum probability of the reaction. The 
cross sections that are governed by Coulomb barriers and resonances show a strong and complex 
dependance on energy and globally decrease steeply towards low energy. 
Cross sections are usually written in the form 

a{E) = ^ e-^-^ (1) 

where S{E) -the astrophysical factor- contains everything that concerns the nuclei and nuclear 
physics and varies slowly with E. The exponential term is related to the Coulomb barrier and 
relative velocity of the nuclei. Experiments which provide measurements of the cross sections in 
the laboratory generally occur at higher energy than the Gamow peak. Extrapolation of S{E) 



to lower E is then necessary. This is difficult to achieve and cannot take into account possible 
unknown resonances occuring at low E. When the S-factor cannot be measured, the reaction 
rates are obtained from pure theory. Recent significant progress in laboratory and theory (hence 
better determination of S-factors down to the Gamow peak) has been achieved as discussed in 
the comprehensive reviews by (23), (24), see also (25). 

2.2.1. Hydrogen burning reaction rates: Uncertainties still exist for the pp chain and CNO 
cycle reaction cross sections in the Sun. They are due to the difficulty to estimate the S-factor at 
low energy and to determine the exact role of the electron screening (both in the laboratory and 
in the star). For most of the concerned reactions, S- factors are determined from extrapolation 
of experimental data to low E but now, for some key reactions, energies corresponding to the 
Gamow peak are accessible through measurements provided by the LUNA experiment at Gran 
Sasso (26). This is the case of the '^He(^He,2p)'^He reaction (see Fig. 1 in (26)) for which 
the S'(£^)-value has not been significantly modified but the error bars are largely decreased. 
Furthermore the cross section of the N^^{p,j)0^^ reaction -the leading reaction in the CNO 
cycle- has now been measured down to energies relevant for IMq stars on the red giant branch 
(see Fig. 19 in (26)) and a resonance is observed. This represents a significant advance although 
an extrapolation is still needed towards solar conditions. Noteworthy, the revision is important 
with a decrease of S{E) by 50%. This has crucial consequences for solar and stellar structure. 
In the solar core, the CNO cycle efficiency is reduced from 1.6 to 0.8% of the total energy. For 
main sequence stars slightly more massive than the Sun the occurrence of a convective core 
depends on the energy production. With the new N^^{p, 7)0^^ rate the convective core appears 
at higher stellar mass and is less massive for a given mass (see for instance the case of a 1.2 
Mq low metallicity star in Fig. 14 by (27)). Finally cross sections that are obtained from pure 
theory can be constrained by helioseismology (28), (29), (23). 

2.2.2. Electron Screening Electron screening is based on the Salpeter's formula (30) with the 
underlying physical picture: the cloud of electrons decrease the repulsive Coulomb effect between 
interacting nuclei with the result of a decrease of the Coulomb potential and an enhancement 
of the reaction rate. This is a static description. It is currently not clear whether dynamic 
effects of the interacting ions can significantly change the impact of the screening in nuclear 
interactions for stars and wherther it must therefore be taken into account. The energy that 
initially fast moving interacting ions have when they get close enough to interact can be lower 
than the mean value of the medium; accordingly their reaction rate is reduced compared to 
Salpeter's prescription. This effect is difficult to quantify and relies on results from numerical 
simulations (31), (32), (33), (34) and the issue is not settled yet. (29) looked at the impact 
of changing the electron screening compared with the classical Salpeter's formula on the solar 
sound speed profile. They found that the solar seismic constraints do not allow variations larger 
than 1% when using GS98. (3) computed a solar model equivalent to the S model but assuming 
the extreme case of no screening at all. Their Fig. 12 shows the difference in the sound speed 
proffie between the seismic sun (12) and model S but switching off e- screening. A decrease 
of the reaction rate by swithcing off the e-screening clearly increases the discrepancy between 
observation and model. This is in agreement with previous results of (23) who showed that the 
large sound speed discrepancy with the use of AGS05 mixture is reduced to the previous level 
obtained with GS98 when the pp-reaction and the screening factor are increased up to 15% 
(Fig. 3 in (23)) In that case, the surface helium and the depth of the convective envelope are 
in agreement with seismic determinations but the sound speed profile in the core significantly 
deviates from the seismic solar sound speed. 



2. 3. Transport of chemicals and angular momentum 

Rotationally induced transport The physical origin of the uniform rotation profile in the radiative 
zone of the Sun unravelled by helioseismology is still debated. Rotationally induced transport 
of angular momentum resulting from a competition between shear-induced turbulence and 
meridional circulation driven by surface angular momentum losses is not able to make the 
rotation uniform in the radiative region of the Sun. For details, the reader is refered to (35). The 
impact of rotationally induced transport on the solar sound speed profile has been investigated 
by (36); (37). (38) have computed a solar model including rotationally induced transport with 
two different assumptions about the initial velocity (slow or 'fast' sun). Their conclusion is that 
for an initially slow rotation, the microscopic diffusion dominates whereas for an initially rapid 
enough rotation, meridional circulation dominates over turbulent shear. However in both cases, 
the discrepancy below the UCZ increases (see Fig. 9 of (38)). This is also illustrated in Fig.l 
which shows the differences in sound speed profile for calibrated solar models computed with 
AGS05 mixture computed assuming either no rotationally induced transport, or rotationally 
induced transport with no surface angular momentum loss or assuming rotationally induced 
transport with surface angular momentum loss. Adding rotationally induced transport makes 
the discrepancy with the observed sound speed worse. Indeed this process smoothes the helium 
gradient below the UZC whereas the observations seem to require an increase of helium in the 
radiative zone below the UCZ. Several prescriptions have been derived for the turbulent transport 
coefficients involved in rotationally induced transport but validitation of these prescriptions 
remain to be done (39). 

Internal wave induced transport in radiative zones must exist in stars as these waves are generated 
at the interface between convective and radiative regions. They have been shown to transport 
angular momentum efficiently enough to make the rotation of the Sun in its radiative part rigid 
(40). One the main uncertainties is related to the absence of a viable quantitative description 
of the generation of waves. 

Magnetic induced transport Instabilities driven by the interaction between rotation and a 
magnetic field could be responsible for the transport of angular momentum and the rigid solar 
rotation (41). This was shown by (37) in the solar case. 

Whether all 3 processes work together to shape the solar rotation profile in the radiative zone 
or only one or two are dominant is not settled yet. 

Investigations of the impact of rotationally and wave induced transport on the structure of 
stars other than the Sun have been performed for instance to explain the Li dip (42); (43). 
(44) studied the impact of including both rotationally induced transport and magnetic field on 
the structure of solar like stars while (45) estimate the seismic consequences and find that with 
the type of dynamo they assume in the radiative zone, the efficiency of rotational mixing in a 
radiative zone is significantly decreased and seismic parameters are then similar to those of a 
non-rotating, non-magnetic star. 

2.4- Near surface layers 

A direct comparison of the observed and numerical frequencies of the Sun shows systematic 
differences that remain small at low frequency but increase with increasing frequency (46). 
Several causes contribute to this offset with more or less importance (47). They are collectively 
referred to as near-surface effects (for reviews, see for instance (48), (4)). 

2.4-.1. Surface turbulent convection One important contribution to the differences between 
the observed and calculated frequencies comes from current modelling of the outer turbulent 
convective layers of the Sun. The description of the convective outer layers of the Sun in ID 
stellar models remains quite approximate due to our inability to represent and implement in a 
ID code a 3D multiscale nonlocal process such as turbulent dynamics (49). At the solar surface, 
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Figure 2. Echcllc diagrams computed with GOLF solar data (blue dots) and using a calibrated solar 
model computed with either an Eddington (left panel) or a Kurucz (middle panel) T(t) law (red dots) 
. Filled symbols correspond to radial modes. Right panel: An additional echelle diagram obtained 
from frequencies of the solar model (red dots) but corrected for near surface effects (green dots) agrees 
with the observations (blue dots) over the fitted interval. 

turbulent convection is inefficient and therefore strongly dependent on the free parameters 
entering the local, ID formulation for convection as well as many other assumptions in the 
formulation. A comparison between frequencies computed with two of the available formulations 
(MLT and CM (50)) for instance shows that the frequency differences increase with frequency 
and reach up to 0.15% at a frequency z/ = 4. mHz i.e. a frequency offset of about 6 — 10/iHz 
((9)). 

Patched models, that is ID stellar models where the outer layers have been replaced in a 
proper way by those obtained with 3D numerical simulations, lead to frequencies in much 
better agreement with the observations (51), (52), (53), (54). For reviews, see (55), (56). 
This approach is valuable for studies of individual stars provided great care is taken in the 
patching procedure, but it cannot be used in a systematic investigation of a large number of 
stars. Indeed 3D numerical simulations with the required quality are quite numerical time 
consuming and therefore not available in the whole range of effective temperature, gravity and 
chemical composition. 

Atmosphere as boundary condition: Another longstanding problem is the description of the 
atmosphere as boundary condition when computing a stellar model (57). In the solar case, 
one can use an empirical atmosphere derived from observations (HSRA) (58) which represents 
accurately enough the Sun atmospheric properties. But in other stars, one must rely on 
atmospheres built using a temperature-optical depth, T{t), law. Commonly used T(r) laws are 
the Eddington law or more realistic Kurucz model atmospheres (59). In order to illustrate the 
impact of the atmosphere boundary condition. Fig. 2 (left panel) shows an echelle diagram built 
from observed solar oscillation frequencies based on SOHO data. An echelle diagram computed 
for a calibrated solar model where the temperature stratification of the atmosphere is assumed 
to follow a Eddington T(t) law is also shown (middle panel). Largest differences are seen at 
high frequency. They are significantly larger than the observational errors. The discrepancy at 
high frequency decreases (roughly by half) when a more realistic Kurucz model is used (middle 
panel) . However one needs extensive grids of such models along an evolutionary track and for 
different masses. Furthermore, model atmospheres suffer from physical imperfections and are not 
applicable over the entire range of needed masses and ages (i.e. gravity, effective temperature) 
(60). 



Nonadiabatic effects: Another source of uncertainty comes from what is usually referred to as 
nonadiabatic effects. They include the effects of interaction of the wave with the radiation and 
with the convection (see for instance (48)). 

All these imperfections concur to generate significant errors on the computed oscillation 
frequencies which in the solar case amount up to 10 /xHz at high frequency. As near surface 
effects cannot yet be reliably included in theoretical frequency computations, an alternative 
has been proposed which consists in removing these effects from observed frequencies. In the 
solar case, this has been quantitatively assessed with the comparison of theoretical and observed 
frequencies of modes with small to large degrees. 

2. 4-. 2. Correcting frequencies for near surface effects In the solar case, due to the large number 
of different modes, near surface effects can be removed from the frequencies. However, other 
procedures must be found for stars other than the Sun. Several approaches have been proposed. 
One consists in comparing the ratio of frequency combinations rather than frequencies themselves 
for low degree modes. In that case indeed near surface effects are at least partially cancelled 
(61). On the other hand, in order to be able to use absolute frequencies, (62) proposed a means 
to correct observed frequencies for near surface effects. They first showed that the systematic 
offset between the observed and theoretically computed frequencies of the Sun is well fitted with 
a power law 

Vm.od{n) = Vobs{n) - a{Uobs{n)/Vmaxf (2) 

with a, h fit to the data and Umax = 3100 Hz is a reference frequency that corresponds to 
the frequency of maximum power in a power spectrum. The frequencies Vobsi^) and Vmodin) 
respectively represent the observed and model frequencies of radial modes with radial order n. 
For the Sun, such a correction leads for instance to small separations (io2 which yield a solar age 
consistent with the meteoritic age (63). 

Fig. 2 (right panel) illustrates the effect of correcting / = frequencies from near surface effects 
according to Eq.2 with r = 1, a = —3.0, h = 4.9 for the same solar model (with Kurucz T(t) law) 
as in the middle panel. The corrected frequency echelle diagram coincides with the observations 
over the frequency interval that was used to derive the parameter values a, h. The next question 
is: how much the parameters a,b, vq do depend on the adopted (Kurucz or another) model 
atmosphere? Of course one must also keep in mind that the values of the fitted parameters are 
valid only over the fitted observed frequency domain (see Fig. 2 left). 

Whether the above procedure can be applied to stars that are different from the Sun is an 
open question. Elements of the answer can be obtained with the use of available 3D numerical 
simulations and ID patched models. Frequency differences for radial modes between a solar 
patched model and a non patched model are displayed in Fig.3 (left) as a function of the radial 
mode frequency of the patched model scaled to the frequency Vmax^o given according to (64) 
scaling law. The curve a{upatch / ^max)^ is also plotted and coincides with the observations over 
the frequency interval that was used to derive the a, b parameter values over an interval up to 
^patch/'^max ~ 1-2. Similar plots are displayed in Fig.3 (right) for 3 other models with different 
effective temperature, gravity or chemical abundance. For these models, we take again Umax 
according to (64) and Tf.ff^Q and Vmax,e from our solar patched model. The parameters a, b are 
adapted to fit the frequency differences over the largest Vpatch/^max interval. Again a power law 
can fit the frequency differences over an interval up to Vpatch/^max ~ 1-2 — 1.4. We note that 
the frequency differences for the 3 models behave differently at high frequency than the solar 
ones. This is likely due to their lower gravity, higher effective temperature. Indeed, for these 
models, the ratio of turbulent to total pressure is higher than for the solar model (the higher 
the effective temperature, the lower the gravity, the higher the ratio Pturb/Ptot) and this results 
in larger differences between patched and non patched models. Note that for the three models 




Figure 3. left: Frequency differences for radial modes Vpatch — Vnopatch between a patched and a non 
patched solar models (blue dots) as a function of i^patch / i^max ■ The black dashed line represents the 
power law according to Eq. 2 with i^max = 3.207mHz,r ~ l,a = —15,6 = 4.9 (For the solar models: 
Teff = 5780, (? = 4.46). right: Same as left for stellar models with different effective temperature, 
gravity as labelled in the left corner. Magenta and red models have a solar metallicity whereas the 
cyan models have a lower metallicity [Fe/H] = —1. Coloured dashed lines represent the power law 
according to Eq. 2 for each corresponding frequency difference. Parameter values arc lymax ~ 1.462mHz, 
r = 1, a = —17, b = 3.9 (red); i^max ~ 1.098mHz, r = 1, a = —10, b = 4.9 (violet); and lymax — 1.071mHz, 
r = 1, a = —15, b = 4.5 (cyan) 

(all three hotter than the solar model), the frequency differences show an oscillating behavior 
in function of the scaled frequency. Whether this behavior is real or artificially introduced by 
the patching process is not known yet. In any case, at lower frequency, the mean variation with 
frequency is well reproduced by a power law up to I'patch/^max 

f» 1.2 - 1.4 

3. Conclusion: 

Despite the huge amount of information provided by helioseismology about the internal structure 
of the Sun, several important open issues remain that we have reviewed. They can also impact 
our understanding and modelling of stars that have a similar structure to that of the Sun. 
An increasing set of such stars are being observed by CoRoT and Kepler, and seismic tools 
developped for studying the Sun are now being adapted to study other stars. 
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